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SUMMARY 


The Surface Target Subgroup of the Joint Technical Coordinating 
Group/Munitions Effectiveness (JTCG/ME) is sponsoring an effort for 
compilation of target test data and vulnerability methodologies. As 
part of that effort, this report has been prepared to provide a ready 
reference to analytical and empirical relationships which are used in 
target \ ulnerability and weapon system effectiveness computer simulation 
routines. The priiicipal intent of this compilation is to provide the target 
vulnerability analyst with an overview of the various methodologies used 
in weapons effectiveness analysis routines which relate'target physical 
and vulnerability descriptions to damage criteria, kill probability esti¬ 
mates, and parametric weapon system lethality analyses. This documept 
should promote a better understanding of how target vulnerability data is 
currently used, and how future target vulnerability definitions may be 
structured to take advantage of improvements in available data and com¬ 
puter technology. The format and level of detail in this report presupposes 
the availability to the reader of the Basic Joint Munitions Effectiveness 
Manual for Air to Surface, Weapon Effectiveness Selection and Require¬ 
ments (Reference 1) and some familiarity, with its contents. 
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SECTION I 
INTRODUCTION 


This report presents a survey of surface target weapon system 
effectiveness computer programs performed for the Target Vulnerability 
Panel of the Joint Technical Coordinating Group on Munitions Effective¬ 
ness'(JTCG/ME). The purpose of the survey was to identify necessary 
input data to the effectiveness programs relating to the target system^, 
so that present and future target vulnerability assessments meet those 
input requirements. The survey showed the existence of numerous 
simulation models, with new models being developed. Most weapon 
system effectiveness models are ad hoc efforts written and used to 
counter a critical emergency or to exploit a crucial opportunity, and 
then filed away to be recalled and modified as required. Most are never 
documented at all; many are documented too poorly to be used by anyone ' 
but the original developer; some, not many, are fully documented. 

Identifying the latter category of well-documented, frequently-used, 
and widely-accepted computer programs for inclusion in the methodology 
survey of this report was a desireable but elusive goal. Many worthy 
programs may have been ignored, preemptorially discarded, or not 
discovered. Many more may have been included in the survey which 
should have been discarded since they offered nothing unique, but were 
modified versions of similar programs or inferior, stripped-down 
replicas of more sophisticated models. 

Of those programs examined, there were basically two levels of 
sophistication: 

• Models containing efficient closed-form or open-end curve 
fitting mathematical expressions, and 

• Models containing inttrative integral solutions such as Monte 
Carlo random sampling + echniques, Gaussian, Hermitian, or 
1-eGendre integration. 

Examples of particular types of programs by name, and recommendations 
for sources of methodology will be made, bnt no attempt will be made to 
catalog all the methodology sources for each type of program since such 
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was not the purpose of this survey. A brief description of the methodology 
available for use by the weapon system analyst and some remarks about 
the degree of appropriateness of such use has been the objective. 

,For air-to-surface target vulnerability and weapon effectiveness 
evaluation, the JTCG/ ME has established standards for computer pro¬ 
grams and has published documentation for acceptable air-to-surface 
computer methodology as a part of the Joint Munitions Effectiveness 
Manuals, Air-to-Surface (Referehcn 2). The methodology includes 
simplified open-end methods, documented in and developed for use with 
the Easic Manual (Reference 1), as well as more sophisticated method¬ 
ology published in a separate methodology volume (Reference 3). In 
addition, a number of well defined special purpose computer programs 
are published as separate documents by JTCG/ ME (for a Listing of these 
programs and a short description of their application, see Reference 4). 
These methodologies are in the process of revision, updating, and im¬ 
provement and are not intended to be the final work in air-to surface 
methodology. The methodologies have been widely used, however, and 
have proved to be of reasonable accuracy for most of the demands and 
applications of weapon systems analysts. 

These JMEM-accepted methodologies are available for use (often 
with many minor modifications) within the weapons effectiveness analysis 
community. They will not fulfill every need of the weapons analyst, but 
they do provide much of the basic framework for developing new computer 
programs for use with unique products of research and development and 
may be considered entirely adequate for establishing weapons lethality 
and force requirements for most inventory munitions. In every case 
where a need fer effectiveness analysis is established, these approved 
methodologies should be considered for use in their pure or modified 
forms and all new superior methodology which must be developed should 
be submitted to the appropriate committees for review and possible 
inclusion within the list of acceptable JMEM methodologies. 

The following sections contain discussions of various methodologies 
for computing target vulnerability and weapon effectiveness. The 
measure of effectiveness or desired outcome of all such computations 
is singly shot probability of damage (SSP/j) or, as it is sometimes called, 
single shot kill probability (P*ss)> which in three-dimensional Cartesian 
coordinates is given by: 




(x, y, z) F(x, y, z) dxdydz 
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where P^r is the total kill probability for ail kill mechanisms from a 
munition detonating at point (x,y,z), and F(x, y, z) is a distribution den¬ 
sity function defining the likelihood of the mumtion detonating within the 
infinitesimal int'. rvals dx, dy, and dz about each point (x,y,z). The 
total kill is a composite of kills from all applicable mechanisms of a 
given warhead combined in the following manner: 





where P%dh is the kill probability from direct hit effects, PkB is the kill 
probability from olast effects, is the kill probability from fragmen¬ 
tation effects, and Pis the kill probability from incendiary effects, 
continuous rod, shaped charge jets, etc. The individual kill probabili¬ 
ties are normally obtained using computerized simulation models which 
are the subject of the succeeding sections. 

The distribution density functions F(x, y, z) in Equation (1) are the 
statistical formulae for measure of dispersion: mean, variance, and 
standard deviation. Most of the programs use a random normal or 
Gaussian distribution for determining individual weapon system ballistic 
dispersion, and since the variance is normally fixed for the 2 (altitude) 
dimension and independent for the other two dimensions (X and Y), a 
bivariate normal distribution is almost universally used for surface 
target effectiveness evaluation programs. For cluster munitions and 
stick releases the distribution may be random uniform witnin the pattern 
intervals or uniformly spaced for intervalomcter releases. Only one of 
the programs reviewed used a modified Kayleigh Distribution for ballis¬ 
tic dispersion, and no other type of dispersion was discovered. 
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SECTION II 


PHYSICAL DESCRIPTION’ MODELS AND DIRECT HIT METHODOLOGY 

.Perhaps the most readily apparent distinction between the highly 
effic.ent simple programs and the more sophisticated and time-consuming 
integral models is the degree of accuracy with which a target element 
or complex can be mathematically simulated for evaluation with the 
several damage mechanisms. As a general rule for the programs 
examined, the most sophisticate^ target models were contained in the 
mor • complicated programs. Older programs were much less concerned 
with realist'c direct hit target models since the direct hit kii> mechanism 
has only recently (with the advent of electro-optical; laser, and infrared 
terminally guided weapons) become a major consideration. A realistic 
physical model in an endgame effectiveness analysis program is only 
considered important to accurately predict damage from a direct hit or 
very near miss of a target element by an individual munition. 

DIRECT HIT EFFECTIVENESS INDICES l • 

Existing vulnerability data for the direct hit kill inechanism are 
quite simple in form. Only three types of direct hit effectiveness in¬ 
dices (El) are currently defined: the vulnerable area (VA), probability 
of damage given a hit (P//#), and bridge effectiveness index (BED. The 
VA of a target is that area which is vulnerable to the damage mechanisms 
of a given weapon. In computing weapon requirements, it is used as a 
mp'sure of the effectiveness for a weapon/target combination requiiing 
a direct hit. VA usually is not associated with any particular segment 
of the target but is a'measure of the vulnerability of thie target as a whole. 
However, the term "vulnerable area" can sometimes be applied to a 
particular segment of «'> target when conditions are suoh that only a partic¬ 
ular segment of the total target is vulnerable to that wjeapon. In that 
case, the probability of damage given a hit (P ffp) is unity over the seg¬ 
ment defined as the vulnerable area (an example is a cfani target whose 
VA would be the floodgate or spillway). 

The number associated with VA is mathematically arrived at by 
dividing the total area within which the weapon can affect the target into 
small units, each with area Aj. Then, the probability of damage given 
a hit on that area (PffDO is found for each unit area, and the sum of the 
products (P// 0 ,){Af) is defined as the VA of the target. 
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Tlu- VA is applied to the calculation of th.' probability of damage 
given a random direct hit on the target (I 1 ///)). The \’np is equal to the 
ratio of the VA to the total presented area of the target. 

For a given munition the Vpp may he a single value which, due to 
lack of data, is usually assumed by the analyst to he either zero or one. 
For those targets for which data exist, a table of P//o values may lx* 
specified as a function of weapon delivery parameters such as terminal 
velocity or elevation angle. 

Tne third form of HI for direct hit kill, the I3E1, is used currently 
for specialized bridge targets only. The BHI expresses the conditional 
damage probability (P///>) as a function of the bridge width (W^) using the 
equation: 

BHI 

- 

PHD = 1 - e 


The BHI is an empirically derived constant which has been developed to 
describe the effectiveness of a particular weapon against a particular 
bridge type. Ppp in this usage represents the statistical probability of 
causing a sufficient level of damage from a random direct hit on the 
bridge roadbed. Similar data for roads, railroads, and aircraft taxiways 
or runways could possibly fill a void which exists in the analysts arsenal 
of effectiveness methodology. 


DIRECT HIT PHYSICAL'SIMULATION MODELS 

Point Location, Circular and Rectangular Models 


There is currently a class of effectiveness models, used not 
infrequently, which totally ignores the possibility of a direct hit on the 
target. In this type of program, the target element or elements are 
defined only by a two- or three-dimensional coordinate point location 
(i. e., the target may be located at (X, Y) or (X, Y, Z), but no size, 
shape, or description is specified for direct hit computations). Examples 
of this type of program are given in References 5 through 13. 
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Only slightly more 'sophisticated.is a two-dimensional model in 
very wide usg within many effectiveness analysis programs. This 
representation of the target is often used for single target elements, 
multielement complexes, or for both elements and complexes within the 
same program. Most of the routines which use a two-dimensional 
physical description model use either a circle or a rectangle in or 
parallel to the ground plane.to represent the target; some pbrmit.a 
choice, and a few even extend the choice to include an ellipse. 

When a circular, rectangular, or elliptical, area is used to repre¬ 
sent a single target element, the radius, (X, Y) dimensions, or major and 
minor axes are normally specified so that the interior of the planar 
figure encloses an area equal to the average presented area of the target 
from the weapon attack aspects. Almost universally, though, the 
practice is to use an estimate of the presented ^rea frorh the 90-degree 
elevation (overhead) aspect without regard to weapon delivery parameters. 
This simplification is necessary partly because tabulated presented areas 
from various'aspects are not readily avaHaole for most target types. 
However, even with exceptional targets for which these data are available, 
the increased accuracy to be achieved in a two-dimensional model hardly 
warrants the effort to include this added precision. References 14 
through 21 are some of the best examples of programs using two- 
dimensional models for unitary targets or complex target areas. 

To represent a complex target composed of many identical elements, 
it is normally assumed that the elements are distributed in a uniformly 
random manner over the entire surface area of the target complex (circle 
ellipse, or rectangle). For this type of simulation, most computer 
models totally ignore direct hit kill of elements. This omission is not 
serious if the area of a target element times the number of elements is 
small compared to the complex area. Even if this is not the case and the 
ignored direct hit kill probability would be large, the error is usually 
compensated for (sometimes overcompensated) by the statistical methods 
for estimating blast or fragmentation kill (see succeeding sections). 

The major attractions of two-dimensional physical models to the 
effectiveness analyst are simplicity, speed, and minimal usage of 
computer time. For simple targets such as a runway, railroad yard, 
warehouse, highway, or bridge roadbed, the model is sufficiently 
accurate. Also, there are occasions when precision is of minor 
importance (o.g., when computing gross trends from large amounts of 
parametric data). The simplicity of the two-dimensional model permits 
very comprehensive simulations to be made in a single computer run. 
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For example, the attack of a target complex consisting of hundreds of 
individually defined circular or rectangular targets by a 13-52 forma¬ 
tion dispersing tens of thousands of submunitions may be simulated 
with such simple models whereas use of a more complicated model 
would be impracticable. The best results for two-dimensional physical 
models are obtained wnen estimating direct hit probabilities for high 
altitude bombing tactics with large numbers of bomblet munitions or 
multiple passes against a single or multipleelement complex. The 
poorest results are obtained when simulating low level bombing tactics, 
artillery delivered munitions, and terminally guided weapons. 

The mathematical formulations for two-dimensional direct hit 
models are obvious for a single target and single munition. For a 
multielement complex target and a uniform distribution of munitions, 
the statistical formula of Equation (4) for F* (fractional kill) is an 
example of the best methodology used in a closed-form program. 



where P/© is the probability of kill given a direct hit on a single target 
element, n is the number of target elements in the complex, Afis the 
presented area of a single target element, Rc is the Radius of the 
circular target complex, and m is the total number of munitions falling 
within the area defining the target complex. 



The next level of accuracy for surface target effectiveness program 
physical description models emplays parallelepipeds or, occasionally, 
spheres and cylinders to simulate the actual configuration of a target 
element. This methodology is almost universal among the high-accuracy 
programs using Gaussian or Monte Carlo integral solutions, and although 
the descriptive accuracy is very poor compared to target vulnerability 
assessment models or air target physical simulation models, it has, in 
the past, been considered by most effectiveness analysts to be wholly 
adequate for surface target weapons lethality assessment of direct hit 
kill capability. Nearly all of the large-scale effectiveness programs 









reviewed which contained a three-dimensional target model use a single 
parallelepiped to model a target element. Each target element is 
described by a point location (vulnerable area centroid), a length, width, 
and height. The most sophisticated large-scale models also define an 
element orientation angle to permit rotation of the target element in the 
ground plane, and permit definition of non-target elements such as 
revetments or nearby structures which would protect or shield the target 
element from a direct hit or fragmentation kill. 

Figure 1 illustrates a typical description of a complex consisting 
of two vulnerable elements and a revetment. Three warhead trajectories 
are shown, one intersecting a surface of each component. The mathe¬ 
matical formulation to determine if an intersection occurs is given 
below. 
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If the defining physical surfaces of a target complex are represented 
as a series of horizontal and vertical rectangles, a direct hit will occur 
if a line defining the warhead trajectory intersects any one of the planar 
surfaces of a target element prior to detonation. This trajectory line 
can be represented by two points (X^, Y w , Z w ) and (Xr. Y r» Zj*). The 
equation of the line is given by: 


X - X„ Y-Y ■ Z-Z„ 
• T T T 



The horizontal or vertical rectangular element can also be 
represented by two sets of coordinates. For example, the first point 
(X/, Y/ , Z/) is located at the bottom end of a vertical surface, and 
the second point (X 2 , Y 2 , Z 2 ) is located at the top of the other end 
(diagonal) of the surface. Using these two points and a third (X j , Yj , 
2.2) derived from them, a plane representing the surface can be formu¬ 
lated in matrix notation as follows: 



Equation (6) defines the vertical plane containing the three points, and 
cap be reduced to the following relation: 






X(Y/Zy + Y 2 z 2 - Y|2^ - Y 2 Z ; ) - YiXjZj + Xfa - XjZ^ 


- X 2 Z ; ) - (X i Y 2 Z 2 + Z/X 2 Y/ - X/Y 2 Z; - ZjX 2 V = 0 (7) 

For simplicity, the following substitutions are employed: 

A — YjZj ~ \ 2 Z 2 ~ Y j Z 2 ~ Y 2 Zi (8) 

B = Xj Z/ + X 2 Z 2 - X/Z 2 - X 2 Z; (9) 

D = X i Y 2 Z 2 + Z;X 2 Y ; - X,Y 2 Zj - Z 2 X 2 Yj (10) 


Thus, the equation of the plane is given by: 

AX - BY - D = 0 (11) 

i 

Solving the equation of the line. Equation (5) for X and Y in terms of Z 
yields: . 

(Z - Zy) (Xjy - Xy) 

X = - + Xy 

Z w - Z T (12) 

(Z - Zy) ( Y w - Yy) 

Y = ---+ Yy 

Z w - Zy (13) 

Substituting Equations (12) and (13) into Equation (11) gives the following 
relation: 






( 14 ) 


a kz^oc*, - x T )(z w ) + imyxz*,- z^a w ) + d(z*- z^i 

A(x ^ _ _ B(Y ^. y r ) “ 


The value of Z obtained from Equation (14) substituted into Equations (12) 
and (13) gives the X and Y coordinates of the intersection of the plane and 
the line. 

The components for the point of intersection are checked to determine 
if they are actually contained -within the bounds of the rectangular surface 
(e.g., is Z; < Z < Z 2 ). Each check is performed by a series of compari¬ 
sons for each coordinate value. 

There are also available a few surface target effectiveness analysis 
programs which more realistically model target elements using a number 
of parallelepipeds for a single target element model. Although these, pro¬ 
grams are currently in use, they are recent developments for evaluating 
terminally guided bombs, and, as such, are almost completely undocument¬ 
ed. Since most are simply modifications of existing programs, it is unlike¬ 
ly that complete documentation will ever be available. The methodology 
for: direct hit determination is identical to that for a single box description 
except that more surfaces must be checked. A typical example of a highly- 
sophisticated three-dimensional direct hit evaluation model is contained 
in the Air Force Armament Laboratory's Program 1721 (see References , 
22, 23, and 24). 


Quadric Surface Models 

Although, among the programs which were reviewed, there were 
no surface target endgame effectiveness analysis routines which use 
quadric surface direct hit models to simulate a target physical confi¬ 
guration, it is nevertheless appropriate to include a brief description 
of this methodology. A surface which can be defined by a second degree 
equation is called a quadric surface. The formulations are available 
(References 25 and 26) and provide much greater accuracy than current 
models although they consume more core space and computer time. 

There is currently much sentiment for and some effort directed toward 
using quadric surfaces in a new family of high accuracy endgame 
effectiveness analysis programs for surface targets. Quadric surfaces 
have long been used for both air target effectiveness evaluation programs 
and for target vulnerability assessment where target physical descrip¬ 
tions are of paramount importance to accurate lethality evaluation. 
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In Lagrangian coordinates, each point (X, Y. Z) on a quadric surface will 
be a solution to an equation of the form: 


AX 2 + BY 2 + CZ 2 + DXY + EXZ + FYZ + GX + HY + PZ + Q = 0 (15) 


Su faces produced by straightline generators (cylinders and cbnes) and 
curves of revolution (spheres, ellipsoids, paraboloids, and hyperboloids) 
are special cases of such quadric surfaces. Other quadric surfaces 
inclv.de regular ellipsoids, hyperboloids, and paraboloids, and combined 
forms such as elliptic cylinders; elliptic cones, elliptic paraboloids, 
and hyperbolic paraboloids. 

All ol‘the above named surfaces have relatively simple closed-form 
defining equations and can be used alone or in combination to define or 
approximate many very complex shapes or figures. For example. 

Figure 2 illustrates how a MIG-21 can be quite adequately represented 
by six ellipsoids. 




Since both the quadric sir face and the weapon trajectory iine can be 
precisely defined in a closed-form mathematical equation, a simultaneous 
solution of toe two equations will define any and all points of intersection 
to accurately predict direct hit of the simulated target component or ele¬ 
ment. An example of the methodology for determining the intersection of 
a line and an ellipsoid follows (from Reference 25). 

Consider an ellipsoid with axes of length 2u, 20, and 27 located with 
the center at the origin of coordinate system B with coordinate axes X*. 

Y b. and Zfl. As shown in Fig re 3, the weapon trajectory line penetrates 
the ellipsoid at points P/ and .* 2 , which are along the X axis'of the weapon 
(P) coordinate system and crosses the Xg, Y b plane at point A. A vector 
from the origin to point A is designated: 



Figure 3. Penetration of an Ellipsoid 





( 16 ) 



20 







z 4 , = sift 9 


(23) 


From Figures 4 and U it can be seen that the coordinates of the points 
P/ and P? are: 


I = - a / x * 


y/ = y# - a / y« 


z i = - a i *« 


where i is either 1 or 2 and a< denotes a scalar quantity, either a. { or a^. 

The general equation of an ellipsoid with reference to the body 
coordinate system is: 


x + y + z 

^ J2 T2 * 1 


where <*, p, and y, are semiaxes of the ellipsoid. Then, substituting the 
Values of x, y, and z from Equations (24) through (26) for the points P* 
and P .2 gives: 

(py) 2 (x B - * t 7L u ) 2 + (»r) 2 (y B - &ty u ) 2 + (*p) 2 (- a,z H )* = (t*py) 2 (28) 


From Equations (21) through (23) and Equation (28), the values of the 
unknown scalar quantity a,- can be determined by solving the quadratic in 
Af and substituting the values of x u , y K , and z„. 


|( B* - 4AC 
2A 


(29) 









where: 


A = cos 2 0 cos 2 * (fti) 2 + cos 2 0 sin 2 * <«'/> 2 + sin 2 0 bfr 

t 

B = 2<ft'f) 2 cos 0 sin * y B -2(/Iy)*cos 0 cos * 

c = (^v) 2 x g 2 + M 2 ya 2 * 


(30) 

(31) 

(32) 


Two values of a . 1. c. a, and a 2 , are found by this method corre¬ 
sponding to points P, and P 2 . If « > 0. then z, < a 2 and: 


B 

2A 


JbTT 


4AC 


2A 


(33T 


a 2 = - 


B 

2A 


i 


B 2 -4AC 


2A 


(34 f 


otherwise, for * < 0. the values of a, and a 2 in Equations <33>_amM34) 


are 


interchanged. The components x„ Y/» and x 2 , y 2 . z 2 can be 
found from Equations (24), (25) and (26). 

Then, from Figures (4) and (5) it can be seen that: 


R/ = X; i + y/ J + z/ 


\ 

k 


(35) 


R 2 = x 2 i + y 2 j *Z2* 

The total length of tue path within the ellipsoid L is: 


(36) 


. Absolute values in Equations (33) and (34) are required to identify the 
entry and exit points in the ellipsoid. 
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4 AC 


L ~ | R, - U, 


A question may be raised at this time about what would happeu .if 
S# : * 11 within the ellipsoid. In this case, the signs of a/ and aj would' 
biff t .mb Equation (24), (25) and (26) would.still remain valid. There 
v.ov.d *he: .r<- he no chance in the solution. 

Idle lution <>t ih.e path of the weapon within the body also involve^, 
by :ef;n;.tu a, a solution lor the angles of impact and exit which the weapon 
ma->*s with the target body [see Figure (6)] . This solution can easily he 
fou-.b_fn.rr. the equation cjf the surface. For a given surfaced, tne grad- 
n-n 7<5>. wil oe a vector perpendicular to the surface at (Xj, y, , z,). For 
an-lhpsoid, it can tne re fore be stated: 




( 38 ) 


where i = 1 for impact and i - 2 for exit. A unit vector N can be formed 
in the direction of the gradient: ' 



Figure 6. impact Angle to Target Surface 
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aru the exit ancle: 

m,'= sin’' k (47) 

G 

wh re: 1 

■> 1 ‘ 7 

(•>•■')*■ > m >» •/-, - (.< ■)* aia ii .sin 4* \’2 1 (■>>)* cos w cos ^ X2 (48) 

V</?V ) 4 x, : + y : : + (oil) 4 zf (49) 

Planar Surface Fitting Models 

Perhaps the most realistic and practical methods for obtaining a 
mathematical description of a physical target configuration for endgame 
direct hit effectiveness analysis programs is the methodology presently 
used in defining target vulnerability data for fragmentation and other ■ 
penetration kill mechanisms. This methodology,, as typified by the 
computer simulation programs SHOTGEN (References 27 and 28), and 
TARGET DESCRIPTION (Reference 29) uses planar surface approxima¬ 
tions for fitting actual physical configurations and can be made as accurate 
as desired by using more or fewer surfaces or points defining the surface 
vertices. 

SHOTGEN uses a series of sequential points such that the third and 
each succeeding point complete the definition of a triangle whose vertices 
are the latt three points in the sequence. By double-pointing or using 
the same point twice in the sequence, a degenerate triangle or straight 
line is defined. The surface of any object, however complex, can be 
approximated in this manner with the degree of accuracy dependent on the 
number of points defined. The more complex the physical shape, the more 
points are required for a given level of accuracy. Figure 7 illustrates a 
target defined by this method. 

TARGET DESCRIPTION uses both quadrilateral or triangular 
planar input to approximate actual physical shapes in much the same 
menner as SHOTGEN. It has no advantage in accuracy since any quad¬ 
rilateral can bo formed from two triangles. The only difference is in 
tin facility for input. The TARGET DESCRIPTION program permits 



26 









input of common figures such as boxes, cylinders, cones, pyramids, 
wedges, etc. , with minimal input. The program automatically breaks 
down the standard figure to quadrilateral or triangular approximations 
and computes, all the vertices. Again, any degree of accuracy can be 
attained if'enough points are used. An example is shown in Figure 8. 

Another program which performs a very adequate function of 
simulating target physical configurations is the COMGEOM (combina¬ 
torial geometry) program which defines the input for the MAGIC target 
vulnerability program (References 30 and 31)., This program uses a 
combination of standard planar figures (parallelepipeds, pyramids, or 
polyhedrons of 4, 5, or 6 faces) and the quadric surface figures, and 
can even define arbitrarily curved surfaces. An example of a tank 
target described in this manner is shown in Figure 9. 

Any of the above three programs can achieve excellent target 
physical simulation for direct hit evaluation and should be considered 
for future development of computer programs for effectiveness analysis 
of terminally guided air-to-surface and surface-to-surface weapon 
systems. 



Figure 8. TARGET DESCRIPTION Model 
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SECTION III 


BLAST DAMAGE YULNEKABI LITY SIMULATION MODELS 


BLAST DAMAGE EFFECTIVENESS INDICES 

To evaluate blast kill mechanisms against various target element 
types or against target complexes, four Els have been defined: 

• Mean Effective Diameter (or Lethal Radius) 

• Crater Diameter 

• Effective Miss Distance (EMD) 

• Blast Mean Area of Effectiveness (MAEft) 

The first El, Mean Effective Diameter, is used to determine the 
effectiveness of mines (see References 32 through 41K It is defined 
as the diameter of a circle, centered at the mine body, within which a 
finite probability of mine actuation and subsequent blast damage to a 
target element exists. This El (interchangeably used or referred to as 
Lethal Radius) is sometimes erroneously used for the evaluation of 
fragmenting mines or active air delivered munitions. In this usage, the 
combined blast and fragmentation kill mechanism effects are used to 
define the Mean Effective Diameter,, or Lethal Radius, and a rough esti¬ 
mate of effectiveness is based on probability of weapon detonation within 
the defined circular boundary. Computer programs using this El are 
statistically optimistic in that the kill probability predictions are greater 
than with other (more accurate) methods. 

Crater diameter is used primarily for linear targets such as roads, 
railroad tracks or airfield runways. In these instances, the usable width 
of the road, track, or runway must be sufficiently cut by the crater or by 
adjacent craters such that passage is denied to traffic. 

The Effective Miss Distance (EMD) of a target is defined as the 
maximum distance from the target edge or surface at which the weapon 
can be detonated and still inflict the desired damage. Calculations which 
consider EMD as the El normally require that the area or volume around 
the target enclosed by the EMD be treated exactly as if it were part of 
the target. , 
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Blast Mean Area of Effectiveness (MAE*) is a casualty or damage 
effectiveness index which is often approximated by a circular pattern. 

The probability of damage given a hit (P hd) within the circular pattern 
radius is unity, and beyond that distance. Pud is aero. The concept 
of MAE & is not used for throe-dimensional simulation. Mathematically, 
the magnitude of MAE $ is computed by the formula: 

r 2 w r R 

MAE* 3 I / Pra <r, e ) rdrdtf (50) 

Jo Jo 

Where P^B (r, 0 ) is the probability of obtaining a kill due to blast damage 
at a distance r and azimuth angle 0 from the weapon to the target centroid 
(r and 8 are measured in the ground plane), and R is a distance beyond 
which Pj (£ is zero for all azimuth angles. Physically, the concept of a 
circular area for representing MAE^ is usually an adequate■ simulation 
technique for evaluation of weapon system blast lethality whenever the 
target is small with respect to MAE (, and impact detonating fuzes are 
used. This is so because blast effects are nearly spherical in their 
volume of influence and because the difference between the distance at 
which 100 percent blast kill is achieved ax d the distance at which no kill 
can be expected is relatively small for a large, range of target types. 


BLAST PHYSICAL SIMULATION MODELS 

Most physical simulation models for blast lethality are relatively 
simple in nature compared to direct hit or fragmentation kill models. 
The simplest of these are, of course, the two-dimensional models. 
There are also more realistic three-dimensional blast models which 
are functionally dependent on distance, and a few in which corrections 
for atmosphere and other dependent variables are made. 


Two-Dimensional Models 


Two-dimensional blast damage simulation models are usually 
extensions of direct hit models, if they co-exist in the same program, 
and are simply concentric similar figures. If the direct hit physical 
dimensions are modeled by a rectangle, the blast model is merely a 
slightly larger rectangle with both the length and width dimensions 
increased by an effective miss distance (EMD) (References 42 through , 
45). 
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If the evaluation model ignores the possibility of direct hit kill, the 
oversight is normally corrected by means of the blast kill formulation. 
The rationale is usually that all direct hits would be well within the lethal 
blast radius and would therefore be scored as a kill by blast. This type 
of two-dimensional blast model is consistent with,the general approach to 
keep evaluation programs simple and efficient; however, for targets 
where direct hit is a critical kill mechanism (bunkers, tanks, armored 
personnel carriers, etc. ), the simple approach does not provide the re¬ 
quired accuracy, especially for terminally guided munitions, low angle 
delivery, or direct fire weapons. If the direction from wliich the weapon 
approaches the target can be considered to be always parallel or perpen¬ 
dicular to the rectangular target area, then for a two-dimensional target 
(i.e., a target without height), the effective target length (L£ 7 -), effective 
target width (Wfp), and effective target area (A^p) are defined as: 

I . . 

I • 

{ 

L Er =L T +2EMD (51) 

W £T = W f + 2EMD (52) 

< 1 • 

*ET = (L ET )m ET ) (53) 

where Lp is the dimension of the target parallel to the flight path and W j 
is the dimension of the target perpendicular to the flight path. For 
example, the target shown in Figure 10 has a target length of 50 feet and a 
target width of 40 feet. If the EMD is 10 feet, then L.£j would be 70 feet, 
W et would be 60 feet and Aft would be 4,200 square feet. 

For a three-dimensional target and an impact detonating fuze, the 
effective two-dimensional target area is a composite of the plan target 
area, the area around the target defined by,the EMD, and the shadow area 
resulting from the projection of the targe? in the ground plane, as shown 
in Figure 11. The shadow length is: 

l SH = tan! (54) 
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Figure 10. Effective Miss Distance About a Rectangular Target 
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Figure 11 . Effective Target Area 





where I is the impact angle and H 7 is the height of the target. The effec¬ 
tive target area is thus: 


k FT = (L r -r 2EMD)(W r + 2EMD) + W r (L SH - EMD) (65) 

To obtain the effective target length (Let) and effective target width (Wet), 
it is sufficient to convert the above irregular projected area into an 
equivalent rectangular area with the width Wet = + 2EMD. Let * s thus 

defined as: 


Ljrr a 


(L r + 2EMD)(W r + 2EMD) + W T (L SH - EMD) 


(56) 


If the target requires a single direct hit to produce the desired 
blast damage EMD is zero. 

For air burst or proximity fuzes it is necessary to include the 
effective target height ( UET ) directly rather than its projection into the 
ground plane (shadow effect), since the lethal blast radius is truly 
three-dimensional. The pseudo three-dimensional technique of including 
the shadow area of the target to determine hit probability of a,three- 
dimensional object in two-dimensions is therefore inaccurate for an air 
blast damage mechanism evaluation model. 

When the MAE* is used as an,El, the effective target length Let 
and effective target width Wet may be defined as: 


L FT s W fr = VMAE b (57) 


The evaluation is always two-dimensional in this event since MAE^ is 
defined as a two-dimensional El. 


Three-Dimensional Models 

The second level of sophistication in blast damage evaluation 
models is the three-dimensional blast volume, within which blast kill is a 
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certainty and 'lor which external' detonation cannot produce a blast kill. 
Nearly one-half of the proa rants which contain three-dimensional blast 
mo els use a simple sphere to model a lethal blast radius around the 
target, and most of tlie remainder use a parallelepiped. Only a very 
few permit the use of Ellipsoids or cylinders,, and none of those reviewed' 
eor'ained more exotic conic sections or quadric surface models. Figure 
12 is an example ol .an ellipsoidal blast envelope about a parallelepiped 
(box) target element. Tins model is a more realistic simulation of blast 
vulnerability than using a larger box to define the least surfaces. 

. Although the target description and vulnerability programs, SHOTGEN 
t’OMGKt'M, and TARGET DESGEI1’TION, are concerned with'describing 
the physical'target and not a lethal blast envelope, logical extension of 
this more sophisticated methodology eould be applied to blast damat 
evaluation. This would require that each of the planar or quadric sur¬ 
faces be projected from a single interior point such that a larger si, ilar 
physical model is defined at the effective, miss distance from the original 
sur'aces. An example of a good, three-dimensional blast model is given 
in AFATL's aircraft shelter vulnerability computer program d'oeumented 
in Reference 46. This program contains an excellent.description of 
blast damage evaluation methodology which is typical of that contained in 
the better three-dimensional blast lethality evaluation models. 
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Blast Damage as a Function of Distance 


The next level of sophistication of blast models is the definition of 
variable blast kill as a function of distance. This method always can be 
reduced to the basic concept of two blast envelopes. The inner envelope 
represents the maximum (usually 100 percent) kill boundary, and the 
outer envelope represents the boundary beyond which blast kUl is 
always zero. Between the two envelopes, some form of mathematical 
interpolation technique defines the probability of bbtaining a blast kill. 
Almost invariably the interpolation techniques are linear, although the 
physics of blast impulse would suggest that interpolation according to 
the inverse square or cube of the distance, or a logarithmic technique 
would be more accurate. 1 


To illustrate the mathematics of such a model, consider a simple 
point target with blast envelopes modeled as spheres. Two radial 
distances are input directly or calculated from input values'cf target • 
critical impulse levels and warhead equivalent bare charge weight of 
explosive. Let Kg/ and Kg2denote these distances as shown in Figure 13 
and interpolation be linear. 



i 


Figure 13. Probability of Blast Kill Versus Distance 


Then the kill probability due to blast, P KB • is g^ven by the following 
relation: 
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1 1 for 0 < D K s/ 

«52' D 

- — for R si < D< *B2 

h B2 ' K fl/ 

\ 0 for 1> > \\ B2 (58) 

where l) is the distance from the warhead to the target. An example 
of this methodology as used in a two-dimensional model for evaluating 
Lethal Area is contained in Reference 14. 

i ’ • , _ . . ' i ) ; 

Atmosphere and Other Variable Dependent Models 

Blast model evaluation methodology in the most sophisticated pro¬ 
grams provides corrections to the size of computed blast envelopes as 
functions of kinematic encounter variables such as atmospheric pressure 
at the target due to altitude above sea level or meteorological changes m 
the environment, or azimuth and elevation angle of the approaching blast 
wave. Since blast impulse propagation is a mass-transfer wave pheno¬ 
mena, the acoustic properties of the transfer medium determine the 
amount of energy transfer. Thus, a given weight of explosive will 
transmit energy better through a denser medium,; therefore, high 
atmospheric pressure and moisture content will support greater target 
damage than a low pressure, dry environment, all other things being 
equal. Since air density varies with altitude above sea level in approv’- 
matelv the ratio ^ for each HI, 000-foot change in altitude, the blast 
impulse I should include a multiplicative correction of: 

' -H 

! = I 5L c~ 31000 (59) 

where Ij^ is the computed impulse at sea level and h is the target height 
(altitude) above sea level. The JMEM/AS open end methods for the Wang 
computer (Reference 17) contains the correction. 

Corrections to lethal blast radii for direction of approach of the 
blast wave are desirable for realistic blast damage probability prediction. 
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Finally, point which is almost universally ignored tor.computing 
bla.-t kill in air-to-surface effectiveness programs is the terminal 
vel city of tlie weapon, at detonation and the effect this has on the forma- 
tun and propagation of the blast wavefront. This effect is perhaps of 
mil m significance for impact detonated 'munitions since the prtmary 
silt ck wave comes f-om a reinforced \ stem effort due to >*e fleet toil of 
the blast wave from, the ground or other surface, For intiueiice and 
proximity fuzed weapons, however, it is of major importance in corn- 
put -ig the magnitude oi‘ an incident blast wave to know if a missile is 
apt 'o.mlung or recoding from the t.ireet at the time .f detonation In 
mc.-t, if'not all, surface'target effectiveness programs, these kinematic 
influences arc not "considered. 








SKCTION IV 


FRAGMENTATION DAMAGE EFFECTIVENESS METHODOLOGY 


FRAGMENTATION DAMAGE EFFECTIVENESS INDICES 

There are only two Els which are used for evaluation of fragmenta¬ 
tion weapon lethality — Mear: Area of Effectiveness and Vulnerable Ar s; 
both have been discussed in previous sections. However, since j 

evaluation of fragmentation kill is often the most difficult and important 
task of the weapons analyst, it is considered desirable to develop these 
concepts further for this section. 


A frank discussion of fragmentation Mean Area of Effectiveness 
(MAE/) should begin with a strong statement to "BEWARE" in the use 
of this El. Unlike MAEj,, there is very little in the way of physical 
realism or tangible simulation concepts which can be attributed to MAE/. 
As a comparison tpol, the concept of MAE is a very utilitarian measure 
of effectiveness. Problems invariably arise, however, whenever MAl|/ 
is converted ,to a more universally understood (or accepted) concept, ! 
single shot probability of damage (SSP^Ifor single shot kill probability 
(Pjtxj)J . MAE/ is a casualty or damage index which is related to the 
average number of target elements damaged by a weapon when the target 
elements (e. g., personnel) are assumed to be uniformly distributed 
over the target area. If targets are uniformly distributed over an area 
with a fixed density of targets per unit area, then MAE/ represents the 
expected number of targets damaged (e.g., expected number of personnel 
casualties). MAKj depends upon the target vulnerability, weapon charac¬ 
teristics. impact velocity, weapon angle of fall, and burst height. For 
evaluating weapon effectiveness against a single target element or other 
than a uniform distribution of target elements, care must be taken to f 
insure that the techniques for averaging effects of MAE/do not negate! 
the applicability of the methodology for the given weapon system. MA|E/ 
should never »- .*ert for SSPp evaluation of directional or focused 
fragment atio > n un; t»ons or conventional cylindrical weapons with low i 
impact angies* against finite target elements. 1 


In the JMEM open end methodology, for the purpose of computation, 
MAE for a weapon/target combination is the area over which a weapon, 
on the average, will cause at least the specified damage to particular 
target elements. Since the effects for bombs and missiles are more ; 
elliptical than circular, am additional factor needed for computations 






the longth-to-width ratio ( a) which can he used to provide a closer 
' approximation to the shape of the damage pattern. To determine effec¬ 
tiveness parameters for this El, an ellipse is obtained whose area is 
.MAE and whose ratio of minor to major axis is (a). The effective 
target length (L£/•), effective target width (W£ 7 -), and effective target 
ar< i (Aet) are therefore defined as: 


1 ET 


= 2 or 1 . 128 VmAE^T) 


(60) 


w fr = Taf 


A fr ~ ^ ,AK / 


(61) 


(62) 


If a is nearly one (i.e., the impact angle is near 90 degrees), then the 
pattern will be nearly circular in shape, and: 



L ET- W ET - 2 


^MAE, 


The mathematical definition of, MAE^ is: 


(63) 


MAE, - f f I*£(r, 6 ) r dr dtf (64) 

Jo Jo 

where P^/r(r, ® ) is the probability of obtaining a kill due to fragmentation 
damage at a distance r and azimuth angle 0 from the weapon to the target 
centroid. A similar formula for MAE*, was given in Section III with 
little or no qualification for its use. The difference in the defming 
functions versus P*£ is important and must be understood to permit 
proper use of the concepts. 

As stated before, for nearly all conventional weapons containing 
high-explosive chemicals, blast effects are nearly spherical (highly 
symmetric) in their realm of influence. The P*- function is therefore 
independent of azimuth angle e. Also, for most targets a sharply defined 
blast damage threshold exists, which means that the P^ function is unity 
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within a threshold radius and drops rapidly to zero outside this radius. 

The fragmentation kill probability function (Pff) is not nearly so 
well behaved. First of all, except for spherical munitions or cylindrical 
munitions delivered at 1Xf elevation, the function is highly 8 dependent. 
Also, there is nothing continuous about a fragmentation pattern composed 
of a finite number of fragments, and since fragments can maintain a let ial 
velocity over a long range, the limit of integration R is a relatively lone 
distance from the detonation point origin. To summarize, the complexity 
of the Pjff function and its discontinuous nature make it unsuited to the 
type* of averaging calculations associated with MAE. The result is a 
tendency to overestimate SSPjj. Given any particular weapon system, 
the standard computational practice of assuming a lethal radius given by: 



(65) 


will normally result in an ovei estimation of SSPp. 

The other fragmentation effectiveness index, vulnerable area (VA), 
is that portion of the target presented area which is vulnerable to an 
individual fragment. The VA is therefore a function of individual frag¬ 
ment parameters (weight, shape, and striking velocity) as well as the 
aspect angles (azimuth and elevation) at which the fragment strikes the 
target. VA is a more basic entity than MAE and, in fact, must he used 
in the computation of MAE. SSP/> for fragmentation can also be cal¬ 
culated directly from VA, and although the integration is more compli¬ 
cated, the results are more accurate. 

Perhaps the best way to define this El is to give a short description 
of the methodology for computing VA using the computer programs 
SHOTGEN or MAGIC. These programs produce the data to compute 
target vulnerable areas as a function of fragment classes (which «re 
determined by a unique striking mass (size and shape) and as a function 
of striking velocity and aspect angles (azimuth and elevation). 

Target vulnerability is described in terms of the vulnerable compo¬ 
nent concept. The substructures ot the target that are most essential 
to the target's continued usefulness, and which cannot be immediately 
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repaired if damaged, are considered to be vulnerable components.. For 
example, if the target is a truck, the engine, transmission,, fuel tanks, 
ignition system, coolant system, electrical system, suspension, and 
possibly the driver may-be considered as vulnerable components. The 
assigned level of damage probability given a hit (P'/®) is specified for 
each vulnerable component and each fragment class. 

The amount of metal that an incoming fragment must penetrate to 
reach a vulnerable component may vary greatly with the weapon's 
direction of approach. For example, a fragment approaching a voltage 
regulator from above might encounter a small fraction of an inch of 
shielding, while a fragment approaching from the front or rear might 
have to penetrate several inches of metal before reaching the regulator. 
The problem is that nothing is known in advance about the orientation 
of the target with respect to the weapon burst point; therefore, a specific 
direction of fragment approach to a vulnerable component cannot be 
specified. At first considerate.n it would seem reasonable to look at the 
components in question from every angle, average the encountered metal 
thicknesses, ^nd consider this as the shielding. The procedure has at 
least one constraint: a weapon may have a penetrating threshold below 
this average thickness. .Therefore, we would calculate a zero probability 
of kill for this weapon even if the weapon can, in reality, penetrate the 
component shielding from certain aspects. 

The SllOTGKN and MAGIC' programs solve the above difficulty by 
computing the shielding thicknesses for each component at prescribed 
azimuth and elevation viewing aspects. The programs set up a grid in 
space which is normal to a ray corresponding to a particular viewing 
aspect. Fragment projectile lines (shotlines), which are normal to and 
randomly located within each cell of the grid, are projected toward the 
target (as shown in Figure 1,4). Each of these lines is examined in turn 
to determine which of the target components is encountered along the 
shotline. The sequence in which component is encountered and the 
equivalent material thicknesses are then stored on a magnetic tape for 
future processing. The smaller the individual grid size, the greater the 
accuracy, and also the greater the computer time required to compute 
all the shotline data for a given viewing aspect. Increased accuracy can 
also be obtained from an increased number of viewing aspects at the 
expense of computer running time. 

After the individual shotline data has been collected, it is processed 
by a vulnerability program such as VAREA (Reference 47)., Besides the 
shotline data, the individual values of fragment classes and initial striking 







velocities must be specified. Each of the ishotlines are examined for each 
targe 1 component in the sequence in which they are encountered and for 
each i agment class and initial striking velocity. This is done to establish 
whether a penetration into a vulnerable component takes place. Pene¬ 
tration is determined with the use of equations (Reference 48) which 
predict residual mass and velocity for fragments impacting specified 
thicknesses of target materials. If the shotline does not intersect a 
vulnerable component or the fragment does not have sufficient mass or 
velocity to penetrate to the component, then no vulnerable area is 
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accumulated for that shotline, For each fragment that does penetrate a 
vulnerable component, the coll area in the original grid is multiplied by 
an appropriate and accumulated as a minute bit of vulnerable area. 

By collecting the vulnerable area accumulated at each variable increment, 
vulnerable area tables are produced which tabulate the data in terms of 
azimuth and elevation angles, and fragment weights and velocities. The 
abeve programs also produce tables of average vulnerable area for indivi¬ 
dual components or for the entire target. This may be accomplished by 
averaging over all possible azimuths at a particular elevation angle. The 
table is then indexed by elevation angle for all elevation angles of concern. 


FRAGMENTATION PHYSICAL SIMULATION MODELS 

All fragmeritdtion effectiveness analysis models which use MAE/ 
as an input effectiveness index are statistical models which operate in 
the same manner as for MAEThese models are efficient closed form 
type programs and are useful for obtaining preliminary estimates, 
processing large amounts of parametric data, or conserving effort and. 
computer time when absolute accuracy is not critical. Examples of two 
of the better methodologies of this type are given in References 4!> and 
50. 


All other fragmentation evaluation programs, including those which 
compute the input values of MAE/ for individual weapons, use VA as an 
effectiveness index and all are basically identical in the methodology for 
computation of P ssK (or SSPfl). An example of this methodology follows. 

Fragmenting Munition Effectiveness Analysis 

Target input data for evaluating fragmenting munitions consists 
of tables of vulnerable areas as functions of fragment striking velocities, 
mass, and aspect angles, as previously defined. The more sophisticated 
models will accept larger arrays representing finer detail in the 
vulnerability model, and they usually will provide for simultaneous eval¬ 
uation of several types of targets and many locations for. each type. 

The fragmentation characteristics of a warhead are defined in terms 
of polar zones bounded by polar angles measured frpm the nose of the 
warhead. Most programs assume symmetry about the longitudinal axis, 
but the more sophisticated programs allow also for variation in frag¬ 
mentation characteristics in longitudinal angular zones as well. War¬ 
head performance is input in tabular form for one or more classes of 
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fragments each having a unique initial velocity, weight, and presented 
area. The tables also define the number of, fragments of each class 
which are ejected from each zone (defined by upper and lower polar 
angles and, if applicable, upper and lower longitudinal angles). 

Other inputs include the terminal lelivery conditions such as 
weapon system approach azimuths, div * angles, deployment altitude, 
number of munitions delivered in a sin ’e pass, delivery errors, ballis¬ 
tic errors, submunition distribution, and submunition terminal velocities. 
The most sophisticated programs have fuzing options for altimeter fuzes, 
active or passive radar fuzes, slant range fuzes, electro-optical fuzes 
and other influence fuzes, as well as impact and delayed detonation fuzes. 
For target complexes containing more than one type of vulnerable, com¬ 
ponent, there must also, be kill criteria specified to determine which 
and how many components must be damaged to constitute a kill of the 
entire target complex. 

Programs using the Monte Carlo procedure select a random 
sample of detonation points for each munition in a single pass and 
evaluate the probability of killing each target component with fragments 
projected from each detonation point. This process is repeated many 
(several hundred) times to obtain average kill probabilities for killing 
each component in a single pass, and the probabilities are then combined 
to obtain the single shot kill probability fer the entire complex according 
to the stated kill criteria. 

Other integration techniques are available which define specific 
evaluation points in terms of standard deviations. For each such point, 
a’distribution weighting factor, is also specified which represents a 
probability that a variable will be found within an interval with boundaries 
between the two adjacent evaluation points. For comparison evaluations 
where only one or two variables are involved in the integration, these- 
fixed tabular integration techniques are generally superior in accuracy 
and efficiency. For three or more variables, the random number 
technique is equal in accuracy and efficiency and considered by many 
to represent a more realistic simulation of the truly random nature of 
weapons delivery processes. 

The mathematical expression defining fragmentation kill probability 
for a single weapon, according t.o the Binomial Distribution Theorem, is: 
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( 66 ) 



where Nf is the number ot fragment classes from the warhead which 
strike the target. Apis the total presented area of the target exposed 
to the fragment hits, Av^ is the area, which is vulnerable to the K f * class 
of fragments, and N* is the number of fragments of class K which strike 
the target. By making logical simplifying assumptions that the distri¬ 
bution of fragments within each class is uniform over a given polar 
(and longitudinal) zone and that the Poisson approximation to the Binomial 
Equation is valid, a simpler form of the basic expression defining P,vf 
is assumed to be: 



- exp 


N f 

• £ p k Av k 
k = i 


(67) 


where P k is the spatial density of fragments from the k^ class striking 
the target (i. e., P k = N k jAp). If the target under consideration is per¬ 
sonnel, the Av* ‘ defined as follows: 



where Apis the target presented area as a function of height of burst 
and horizontal distance, and is the kill probability based on a 
fragment hit of a certain mass and velocity. The presented areas for 
standing, prone, and foxholed troops in various types of terrain as 
defined by the Ballistic Research Laboratories (BRL) may be input or 
compiled into the programs. The BRL conditional kill probability 
expression (Reference 51) is used to compute p cp 

To evaluate Equation (67) the following information must be deter¬ 
mined: 
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• Each class of fragments from the warhead which strike the 
target 

• Striking velocity of the fragments of each class 

• Striking aspect angles of the fragments of each class 

0 Striking density of the fragments of each class. , 

Some programs assume a static situation (i.e., that the warhead is on 
the ground and stationary at the time of detonation), but the more accurate 
programs calculate relative velocities and dynamic fragment character¬ 
istics. The dynamic fragmentation characteristics are obtained by- 
vectoring the static warhead fragmentation data forward with respect to 
the missile velocity vector. Hence, each static class of each polar zone 
is redefined to form a new dynamic fragment class with its own upper 
and lower polar angle limits. 

These dynamic polar angle limits for each class of fragments are 
used in conjunction with the following relations to determine which of 
these classes have fragments that strike the target at the point location 
<x r , y-f, 7-j) under consideration. The angle from which a fragment 
must emerge to strike the target is given by: 



where 6 is the angle formed by the missile longitudinal axis and a line 
between the detonation point and the target centroid, (x w , y w , z w ) is the 
detonation point coordinate of the warhead a and g are the warhead 
terminal flight path angles at time of detonation and D is the distance 
from the warhead to the target given by: 

D = £(x w -x r ) 2 + (y*,- y r ) 2 + (i w -7. T ) 2 ~^ i (70) 

The resulting fragment striking elevation angle is given by: 

? = sin -/ £(z w -Zt-)/dJ (71) 
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At this point the distance to the target and the fragment classes 
which Strike the target are known* T. calculate the atriWogdena^ 
of fragments of each class, the initial dynamic velocity <V 0 > must 
calculated as follows: 


V 0 - V m cos o 1 (V s , 


? 2 
„ sin $) 


where Vm is the missile terminal velocity. V„ is the sUUcit^Ui velocity 
for the fragments of the k<* class, and ♦ is as computed m Equation (6-)., 
Finally, the dynamic fragment striking density, P/y is compu. 


n 2 v,. 




where > H is the static fragment density (fragments per storadian) of the 
k lh class* and \' 0 is defined by liquation (*2). 

The remaining tasks are to solve hy the 

the* kill probabiiity due to fragments 
(Pgf), as defined by Equation (67) may now be calcula*ed. 

It should be reiterated at this point that 
vulnerable -e a concept^ component is 

for the VA to be partitioned among several point locations rath 
a single one. 

Most of the effort to develop better concepts is being <^cted 
toward improved air target effectiveness programs where > the> ext ent 
of ‘he concept inadequacy and hence the need for improvement g 
Howaverone new concept for development of superior surface forge. 
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effectiveness methodology is the extension of the POINT BLKST 
vulnerability assessment routines (e.g., Reference 52) to the calcula¬ 
tion of endgame effectiveness of surface targets. This effort is 
currently being conducted at AFATL. 


Multiple Fragment and Synergistic Effects 

Fragmentation kill probability (P xph as.computed by the above 
methodology, is restricted to the precept that fragmentation damage is 
invariably the result of a single fragment penetration into a single 
vulnerable component. There is no allowance for the real-life possi¬ 
bility that one fragment might penetrate the shielding and a second 
(following) fragment could damage the vulnerable component behind the 
shield, or that the simultaneous impact of many fragments might cause 
damage in excess of that which would be caused by the same number of 
fragment hits, one at a time. Only recently has much thought been 
given to the possible bo. is effects of multiple fragment or synergistic 
damage mechanisms. Recent deve^pments in the focusing of kinetic 
energy fragments, primarily for antiaircraft weapdns, has added an 
impetus to u '’ Mop realistic methodology to account for these recognized 
phenomena. Thus far, most of the method viogy. has concentrated on 
redefining the vulnerability to include essential structural components 
which were previously considered only as shielding. The damage 
mechanisms are basically of two types, "Energy Density" or "Material 
Removal. " For the former, a threshold energy density is defined and 
if the kinetic energy of all the fragments within a unit area collectively 
exceeds this value at the time of impact the target is sufficiently 
damaged. This concept assumes that all of the fragment kinetic energy 
is transferred to the target (the fragments do not penetrate). The 
second concept is almost the antithesis of the first in that it assumes 
complete penetration of each fragment, removing the structural material 
from each fragment hole and infers that the damage results from weak¬ 
ening the structure sufficiently so that normal static and dynamic loads 
cause complete failure. Although much attention has been devoted to 
these concepts, the resulting methodologies are not yet widely accepted 
and are seldom, if ever, used in evaluation of surface target weapon 
lethality. 





SECTION V 

MISCELLANEOUS SPEC IAL PURPOSE VULNERABILITY 

methoixhjOGY 


The previous sections on direct hit, blast, and fragmentation 
vulnerability were attempts to describe the methodology available in 
mere than 90 percent of the endgame effectiveness analysis computer 
programs. There are perhaps three other conventional kill mecha¬ 
nisms worthy of a brief mention. These include incendiary/fragmentation 
weapons, shaped charge weapons, and flechette weapons. The following 
paragraphs present a very cursory discussion of the most significant 
methodologies which are unique to these particular' 1.11 mechanisms. 
There are also several continuous rod warhead types in the inventory, 
but these warheads are used exclusively for air targets and will not be 
included in this surface target report. 

INCENDIARY/FRAGMENTATION MUNITION MODELS (REFERENCES 
53 THROUGH 58) 

The most inclusive treatise on effectiveness of incendiary or 
incendiary/fragmentation munitions examined is contained in the com¬ 
puter program documentation of Reference 56 and 57 which presents a 
description of a computerized spatial-time model of events potentially 
leading to ignition of fuel vapors by exposure to a mix of fragmentation 
and incendiary munitions or to munitions containing burning reactive 
metal particles. Thus, vehicular targets (e. g., trucks) containing 
fuel, become Considerably more vulnerable when incendiary materials 
or reactive metal particles are added to the conventional fragmentation 
munitions. 

The salient events occurring during an encounter between a reactive 
fragmentation munition and a target with flammable fuel are as follows: 

• The munition or munitions detonate and project large primary 
fragments and numerous small burning particles of incendiary 
or reactive material. 

• The primary fragment perforates the fuel container wall, 
passes through the fluid, and causes a cavity to expand in 
the fluid. 






• The cavity begins to collapse half-way through its period and 
forces a highly atomized spray of fuel through the perforation 
during completion of ns period. 

• The fuel spray travels in a conical pattern away from the 
container wall as tin* cavity repeats successive pulses whose 
volume and duration are degraded until a steady leak results. 

• j The burning particles can internet with the fuel spray or with 

the resultant puddle of leaked fuel and cause ignition if the 
interaction occurs after spray emergence and prior to 
particle burnout. 

J ' 

• | The ignition can propagate to the main fuel discharge and a 

sustaining fuel fire results. 


Since gasoline is highly flammable at most ambient temperatures, a 
sustained fire is readily achieved due to the surplus of vapors in the 
local ared. However, diesel fuel does not vaporize at most ambient 
temperatures and, although igniting the highly atomized fuel spray by 
burning particles is easily accomplished, a sustained fire is much more 
difficult ta achieve. A considerably narrower lethal time span for sus¬ 
tained ignition of diesel fuel by burning particles exists than for Mogas. 

In particular, it appears that a significant delay in ignition time enhances 
probability of sustaining a lire to the main discharge stream since the 
spray is fully exposed, :t delivers a larger tnormal pulse, and the spray 
is not consumed prior to main stream issue. 


Burning ['article Protection 


Static particle projection charaetoristics are (as with normal frag¬ 
ments) defined in polar zones, each of which may have distinct values 
for: ! 


o I Upper and lower polar angles limiting distribution (♦yand ^j). 

i ‘ 

• j Number of particles (N p ). 

• Lethal burn time (t^). 

• I Initial velocity (V#). 





Drag coefficient limits (K/ and Kd. 


The burning particles are assumed to be launched in spherical segments 
formed by and Qq. They obtain a constantly increasing radial iis^ribu- 
ticn during flight due to uniform distribution of drag factors between 
K / ind K? . Therefore, at any time from detonation, he burning parti- 
i'l< s are located in a solid spherical segment defined by the polar angles 
and outer and inner radii: 


R 0 (t) = 1/K, In (K/ v, t + 1) , (74) 

R; (t) = 1/Ky In (K? v 0 t c 1) (75) 

T’u d Container Perforation 

The large primary fragments launched at munition detonation, 
impact and perforate the fuel container. The Project THOR equations 
(Reference 43) are used to determine perforation capability and resiauql 
mass and velocity of these fragments. The maximum azimuth deviation 
from the container surface normal is related to impact velocity (v), 
mass (m), elevation angle (ELKY),.and wall thickness (T,) as follows: • 

i 

AZ max = AZ U, m,T f> ELEV) (76) 

The residual velocity (v^) and mass .(m/p) were also related to impact 
characteristics as follows: 


V R 1 

v* <v, ni, $, 

T„) 

dps) 

(77) 

m R 

m R (v, ra, 9 . 

T e ) 

(gr) 

(78) 


"The specific functions arc not identified in this text to avoid security 
classification. For further detail see Reference 57. 
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where 0 is obliquity, given b\: 


6 cos ^ £'*os (AZ) .'os (KLEV)} (79) 

Tho fragment is assumi'd to cause a perforation size equal <.s pre- 
si'ntt'd area (A p). 

U ~ uel Cavuation and Spray Knmrprncc 

After the primary fragment has. perforated the fuel container wall, 
tlie residual fragment begins to pass through the fluid. A cavity is 
formed in the fuel behind the residual fragment. The volume and periou 
of the cavity has been shown to vary with fragment kinetic energy as 
follows: 

V f = 8. 92 x 10' 7 E (cc) (80) 

T s <l. 8f) x 10' 6 t}R (sec) (31) 


where E is fragment kinetic energy in ergs obtained from: 

E^Jim^v^- (82) 


Equations (80) and (81) were derived from underwater ballistic experi¬ 
ments. 

Assuming the cavity achieves its maximum extension at one half 
of its period, then the fuel spray begins to emerge from the fragment 
perforation at time r/2 seconds after impact when the cavity begins to 
collapse. Emergence velocity of the spray, assuming uniformity, is 
given by; 

v, = 0.30 V f /[{r/2) <A p )] 

r>3 


(83) 







Th ' spray emerges from the container in a conical pattern normal to 
the surface within a cone angle of approximately 20-degrees. The outer 
and inner diameter limits of the spray Dq and 0/ are expressed as a 
function of time (t) from detonation as follows: 



D () (t) 


Oll t v! t * r/2 

(I/k) in [Kv, (t - t. - r/2) 4 l] 


D f (t) = 


0 if t < tj- + r 

<1 /K> In [Kv, (t - t f - r ) + l] 


(84) 


(85) 


where K is drag eoeffieienl of the fuel droplets, and tf is the primary 
fragment time of flight, given In : 


t ; = (V KfR - 1 ) / ( Vy K f ) 


( 86 ) 


where K/ is drag coefficient for the primary fragment, R is distance 
from detonation point to the fuel container, and vy is initial velocity of 
the primary fragment. Equations (84) and (85) define location of the 
first fuel spray in front of the container in terms of time from detonation. 

It should be noted that wall thickness (T e ) of the fuel container is 
the only target parameter in these equations. Fuel spray characteristics 
are sensitive to other target parameters including container size, 
container shape, ullage/fuel ratio, impact point relative to fuel head 
as well as fragment configuration and mechanical properties. These 
relations are based on full 5-gallon military fuel cans impacted at the 
center by steel fragments. 


Fuel Cavitation and Spray Pulses 

The initial fuel cavity collapses and begins to expand again with a 
somewhat degraded maximum volume and period. As many as five 
consecutive pulses of fuel spray have been observed during tests against 
full 5-gallon fuel cans. However, three pulses were observed m most 
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cases; thus, the model enables the program user to evaluate effects of up 
to three pulses. 

The time from impact when the first fuel pulse issues from the 
container is given by: 

. . t, = Tjl 2 (87) 

where tj is the period of the first pulse, given by Equation (81). 

The empirical approac h nvolved establishing a value for a which 
forces the following equations.'to fit observed issue times for the second 
and third pulses: 


t 2 = Tj + " Tj / 2 

(88) 

tj - Tj 4 « T J + <T- T;/2 

(89) 


This establishes a degradation factor for the period of each consecutive 
pulse. Additionally, the cavity volume at maximum extension was 
degraded as follows: 

V 0 = V«.,/(l + «/2) (90) 

Z i 

\' Cj = V (1 + + «72) , (91) 

where V CJ is maximum volume of the first cavity [Equation (80)], and 
\J., and V C j are volumes of the second and third cavities. A value for 
«»of 0.7 has been derived from high speed film of gun tests. This value 
provides good correlation with observed times of 1 2 and tj, and also 
correlates well with spray issue velocity calculated by Equation (83). A 
space-time history of the second and third pulses is established by 
recycling through Equations (83), (84), and (85) using modified data for 
cavity volume and period. 
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If the primary projectile perforates the container beneath the fuel 
level but near a free surface, the cavitation phenomenon is not expected 
to behave with no edge effects as described above. However, the nature 
of events involved in the cavitation process are of such complexity as to 
preclude their formulation from a purely theoretical standpoint. 


SHAPED CHARGE MUNITION MODELS 

A second class of special purpose surface target programs consists 
of shaped charge simulation models for both conical and linear shaped ' 
charge munitions. The shaped charge munitions take advantage of the 
focusing properties exhibited by a detonation wave encountering a cavity 
or channel which is normal to and symmetric with the advancing detona¬ 
tion wave front. Figure 15 illustrates both a conical shaped charge and 
a linear shaped charge. The primary kiil mechanism for both weapons 
is an ultra high-velocity jet of molten metal which is highly efficien* in 
penetrating structural materials or armor. 1 ■ 

i 

A direct hit by a conical shaped charge munition will produce 
penetration of from 4 to 6 calibers through hardened armor plate. This 
means that a 6-inch-diameter munition will consistently penetrate from 




Figure 15. Shaped Charge Geometry 
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24 to 36 inches of armor. Most direct hit programs for conical shaped 
-charge munitions assume a penetration probability of unity and merely 
compute the probability that a randomly located penetration shotline 
will impact a vulnerable component. Reference 59 is an example of a 
highly accurate conical shaped charge evaluation program. 

Linear shaped charge evaluation methodology is more complicated 
since there ai*e several separate kill mechanism sources, and each is 
a line source of high velocity, molten-metal jet rather than a point 
source. An example of a typical formulation for a linear shaped charge , 
evaluation program follows (from Reference 60). 

In defining linear shaped charge warheads for effectiveness compu¬ 
tations, it is necessary to evaluate three kill mechanisms: blast, jet- 
slug, and fragmentation effects. Each kill mechanism is capable of 
inflicting kills on targets by severe structural damage. The target 
structure may be represented as previously discussed in Section 11. 

A linear shaped charge warhead will have blast kill potential 
equivalent to a slightly smaller conventional warhead except in the 
jet-slug zones, where the blast kill will be enhanced by increased 
overpressure and consequently higner blast kill capability. To model 
the shaped charge blast kill, the blast ellipsoids are defined and checked 
initially for conventional blast kill. If the target is not killed by con¬ 
ventional blast, a second set of larger blast ellipsoids is checked to 
determine if the target is within lethal range of the enhanced blast. The 
final check will be to determine if the blast center is contained within 
the warhead jet-slug zones. The probability that the target is killed by 
the damage mechanism is called P g. 

A linear shaped charge warhead also projects slugs or rods within 
the jet-slug zone oriented parallel to the warhead centerline. Up to 
• three of these slugs for each jet may be projected. The slugs are 
described by mass, cross section, ballistic data, initial velocity and 
flyoff direction. The slugs are projected along the centerlines of the 
angular segments representing the jet-slug zone. Target vulnerability 
to slug impact is defined by a series of ellipsoids representing structural 
components. Encounter geometry, slug velocity, drag, and direction is 
used to establish the event of a slug-target intercept with the target 
ellipsoids. Each ellipsoid has a conditional kill probability assigned 
for si*' \ impacts. The probabilities are combined statistically, assum¬ 
ing independent events, for more than one impact into the probability of 
kill due to jet-slug (Pj). 
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In addition to slugs, fragments with a spectrum of masses, initial 
velocities, and drag characteristics are projected from each linear 
shaped charge liner. The capability for representing several classes of 
fragments is provided. The fragment spray is limited by two angles 
measured from the forward and aft ends of the warhead, and an angular 
increment about the warhead centerline. 

Targdt vulnerability is represented by a vulnerable area centered 
about one or more points. If a vulnerable point is established to be in 
one of the fragment sprays, then kill probability due to fragments is 
computed by: 

N 

P = 1 - EXP < -Xj f>i AV,) (92) 

whore p,- is impacting fragment density of the i tfl class, AVj is target 
vulnerable area of the i fA class, and N is the number of fragment classes 
striking the target. The computer program formulation includes the 
capability to input tabular vulnerable area as a function of fragment mass, 
velocity, and striking aspect. The capability for computing structural 
kili from multi-fragment impact is also included. 

The basic geometry for a specific weapon-target encounter, includ¬ 
ing aim point, guidance error, fuzing point, and detonation point, is 
established identically for linear shaped charge warheads as for con¬ 
ventional fragmentation munitions. However, techniques for computing 
fragmentation kill probability for conventional warheads, assuming a 
point source of fragments, are not adequate for evaluating linear 
shaped charge warheads because of the narrow polar beam spray zone 
(on the order of 5 degrees) which includes all fragments projected by the 
warhead. For example, the width of a 5 degree beam spray at a distance 
of 50 feet from the point warhead is < 5 feet, thus a warhead length of 1 
foot, would constitute an error of 20 percent in the actual beam spray 
width if a point source of fragments were assumed. The following para¬ 
graphs present techniques for correcting beam spray zone width to 
account for warhead length. 

/ 

The fragmentation subroutine for the conventional warhead computer 
program starts to compute fragment kill with a target description and a 
warhead detonation point (X^, Y D , Z D ). The point of fragment-target 
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intercept [X, (i, j), Y, (i, j), Z, (i, j)] is found which determines the 
points where fragments projected from the warhead with velocity i hit 
the target vulnerable component j. For.an intercept, the polar angle s^ 
from which the fragments of the warhead were projected to strike the 
target component is computed. This angle is then compared to the for¬ 
ward and aft beam spray angles Bj and 0 2 to determine if the solution is 
valid (0/< 0f <02)- Fbr a linear shaped charge, the values of 0/ and 62 
must be modified tp account for warhead length. 

Figure 15a is a diagram of the current fragmentation model for a 
cylindrical warhead showing the warhead detonation point (X D , Y D , Z^), 
the point where the fragments strike.the target vulnerable point (X*. Y,, 
Zj), the forward and aft polar angles 0; and 02, and the angle to the 
target Bp. V. .th this model it is obvious that the vulnerable point is out¬ 
side the mathematically defined beam spray zone (0p< 0/) and the solution 
would be rejected. Figure 16b shows how the target vulnerable point cam 
be included in the beam spray and establishes the geometry for formu¬ 
lating equations to correctly model the fragment-target intercept. 

The formulation is as follows: 

From Figure 16a the distamce of fragment flight Df is given by: 

a 

D F = [(X, - X^ + (y, - y D ) 2 + (z, - z D ) 2 ] (93) 


The perpendicular distance L p from the warhead centerline to the target 
vulnerable point is given by: 


L p = (D F )(sin Bp) (94) 


The width of the forward beam spray at the perpendicular distance L p 
is defined in the point source model (Figure 16a) by: 




tan Bj 


(95) 
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(a) Point Source (Uncorrected) 



Figure 16. Fragment Spray Representation 


6 


s 




The width of the aft beam spray at the perpendicular distance Lpis de¬ 
fined by: 



S 


2 ~ 



tan 8 j 


(96) 


Adding one-half the warhead length (W^) as shown in Figure 16b to each 
value gives: 



(97) 


S *= - 


W, 


tan 0 2 2 


(98) 


The corrected angles Oj and 0 2 are then given by: 


8 j'- tali' 1 


hs 

S /' 


(99) 




( 100 ) 


If Op is'between the corrected values 0j ' and (0/ *< Of < 0 2 ) the vulnerable 
poinf will be in the fragment beam spray. 


The density of fragments which strike the target must also be 
established to compute kill probability. Hence, techniques employed 
for point-source fragmentation cannot be directly applied for evaluating 
a linear shaped charge warhead. The equation which determines the 
dynamic striking fragment density for conventional fragmentation 
warheads can be employed with two changes: 
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• Use the modified polar angles, dj‘, and to account for 
warhead length. 

• Multiply by a factor of 2 *7AV, where At is the roll angular 
increment containing the fragments projected from a liner, 
to reflect the fact that a linear shaped charge warhead is not 
isotropic. 

Kill probabilities due to fragmentation (P f) are determined account 
ing for fragment drag, impacting density, mass, velocity, angle, and 
dynamic missile velocity. 1 

Thus, all three kill mechanisms for a linear shaped charge are 
accurately defined and the total kill probability (P) resulting from blast, 
jet slugs, and fragmentation is then formulated as: 


p = i - (i - p b ) a - p s) a - pf> aoi) 


FLECKETTE PROJECTING MUNI 
61) 


ION MODELS (REFERENCES 60 and 


Warheads which project flechettes as a kill mechanism may be 
assumed to be special cases of fragmentation warheads. Flechettes may 
be thought of as stabilized, streamlined fragments which travel through 
the airstream like darts or arrows presenting little aerodynamic resis¬ 
tance and thus retaining their velocity over long distances. 

, i 

The major differences in evaluating flechettes and fragments is in 
the distribution patterns and trajectory ballistics. Flechette patterns 
are severely limited in the maximum divergence angles, and because 
they maintain a lethal velocity for (long distances, the effects of gravity 
cannot be neglected. The following paragraphs present an example of 
flechette methodology for personnel target lethal area computation 
from Reference 59. Similar program formulations are available for 
material and mixed 'personnel-material) targets. 


The warhead is assumed to release flechettes in a conical pattern 
with prescribed angular limits defined by a beam spray half angle 
measured off the longitudinal axis: of the warhead. The flechettes can be 
distributed within this conical pattern in a uniform or normal manner, 
determined by an indicator in the (input data. Beam spray half angles a re 


i 


52 






limited to 30 degrees or less and the warhead is assumed to be sym¬ 
metric about a vertical plane containing the longitudinal axis. The 
conical pattern defined at release is subdivided into increments along 
the beam spray half angle, measured off the longitudinal axis of the 
warhead, and further divided into radial increments measured around 
and perpendicular to this axis (Figure 17). These angular increments 
of the flechette pattern are then projected into the ground plane by 
trajectory simulations. The program is capable of simulating a war¬ 
head containing up to three distinct classes of flechettes defined by: 

• Flechette mass. 

• Initial flechette velocity. 

• Beam spray angle and angular increments. 

• Distribution and number of flechettes. 

• Flechette reference diameter. 

• Cp versus Mach number curve. 


In addition, up to five sets of flechette reference diameters and/or 
C/) versus Mach number curves may be used as a function of real time 



Figure 17. Sample Warhead Geometry 
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of flight for each class of flechettes. This capability enables the simu¬ 
lation of unstable fleehette flight. 

The Fleehette Lethal Area Program was compiled to evaluate the 
lethal area of fleehette projecting warheads against personnel targets. 
Munition lethal area is established by integrating kill probability against 
either standing or prone troops'or troops crouching in foxholes within 
the impact pattern area. Vulnerability of personnel targets is defined 
in terms of target presented area and a conditional kill probability 
established for a ,hit on the target. Presented area of the target is a 
fut ction of the orientation of the fleehette at impact. Conditional kill 
probability is established through a casualty criterion equation which 
utilizes fleehette mass, striking velocity, and three constants defining 
the kill criterion (see Reference 61). 

Terminal delivery conditions of the missile are significant in that 
a damage assessment point in the lethal area grid must be determined 
by integrating the equations of motion from specified initial conditions. 
Velocity and direction of a given fleehette at release predetermines the 
point of impact in the ground plane. The warhead, as discussed above, 
is subdivided into a grid defining orientation of a fleehette dispersion 
pattern, where each point in this pattern is integrated through time to 
establish an impact pattern in the ground plane. Release geometry at 
launch is defined relative to the missile coordinate system. 

Initial velocity (VM) imposed on flechettes is determined relative 
to the ground plane. This is the total velocity of the fleehette and 
represents the vector sum of the missile velocity and any ejection 
velocity associated with the dispersion mechanism. At the time of 
release, stability of the fleehette becomes a problem that is not readily 
defined. However, the tumbling of an unstable fleehette may be simu¬ 
lated by adjusting a portion of the drag curve where instability would 
occur. Since the reference diameter is directly proportional to the drag 
it may also be altered to account for instability. Up to five sets of 
fleehette diameters and/or Cp versus Mach number curves may be input 
for each class. These data are a function of time of flight and enable 
the simulation of turn., a " t , oscillating, and stable motion. A smaller 
time ' wrement of integration may be employed in the trajectory compu¬ 
tation for a specified period of tine' following a change from one set of 
drag data to another. Th’ , is to ensure accuracy in the iterative solu¬ 
tion of the equations of motion. 

Height of release is defined by input data and is common for all 
classes of flechettes. This program was especially designed to cope with 
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high altitude release, as the gravity effects are more influential for 
this terminal release condition. 

The warhead terminal attack angle (measured positively downward 
from the horizontal to the nose of the missile) is input to relate the 
flechette flight direction to the ground. This is necessary as the dis¬ 
persion pattern is defined relative to the missile axis, and some relation¬ 
ship must be established between its orientation and the impact plane. 
Gravitational forces are assumed to be normal to the ground plane. 

The methodology contained in the Flechette Lethal Area Program 
was formulated from basic mathematical and physical equations. It may 
be subdivided into three orderly sections discussed in detail below. 


Warhead Methodology 

The simulation of a flechette dispersing warhead is formulated 
by geometrically establishing an ejection grid, relative to the warhead 
or missile axis as discussed above. Each point in this pattern isj 
associated witu the directional vector of a flechette at the time of release. 
The vector is assigned an initial velocity obtained from the input Idata 
and is transformed into its three components along the coordinate axes 
of the ground reference system. These three velocity components given 
by VX, VY, VZ are then used in solving the equations of motion for the 
trajectory routine. 


Trajectory Analysis (Figure 18 ) 

All effectiveness methodology currently available contains the 
assumption that fragments fly straight line trajectories from warhead to 
the target. This assumption is valid if velocity is high, distance) to the 
target is low, and hence time of flight is small. However, flechettes 
launched at high altitudes have a relatively long flight time; thus; gravity 
significantly affects the trajectories of the flechettes. To accurately 
account for these effects the equations of motion must be integrated from 
specified initial conditions (velocity and direction) at launch, through 
time and space, to impact. The lethal area integration grid must be 
established by initial conditions of projectiles at the warhead point of 
detonation. The equations of motion defining the flight of each flechette 
are derived from Newton's second law where the drag force is given by: 
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Y d s '1/2 C d p AV ■= - 


(102) 


VtLOC/nr OF MISSIll 



Figure 18. Trajectory Geometry 


Expanding this formula into component forms relative to the ground 
reference system, the equations of motion for a given flcchette can 
generally be described by: 


X = - (1/2 C D P A/M) (VHV X ) 

Y = - (1/2 C D P A/M) (V)(V,) 

Z = - (1/2 C D P A/M) (V)(V r )-g 


(103) 

(104) 

(105) 


where C^)is the drag coefficient, P is the density of air, A is the pre ¬ 
sented area of the flechette, V x , X’y, V x are the velocity components of 
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the fleehette velocity vector V, M is the fleehette mass, and g is the 
acceleration of gravity. 


Lethal Area Computations 

Integration of Equations (103), (104), and (105) can. ot be performed 
in closed form and requires a finite difference iterative procedure. The 
numerical solution of these equations is performed in the program by the 
method of Rungc and Kutta. From the solution of these basic equations of 
motion, the striking velocity, impact angle, and the impact coordinates 
(relative to the point of release) of the flechettes in th* ground plane can 
be determined. The impact pattern is determined by looping the above 
procedure for each point of release in the ejection grid. The polygons 
formed in the ground pattern by the impact coordinates are sequentially 
numbered as shown in Figure 19. 



OUTPUT POINT NUMBERING SEQUENCE (J.K) POLYGON NUMBER 
SEQUENCE. N; WHEN N = 1.2.n 


Figure 19. Sample Impact Pattern Geometry 


As discussed above, methodology has been developed to evaluate 
munitions projecting up to three classes of flechettes. This makes it 
possible for the simulation to contain one, two, or three impact patterns 
which may or may not overlap each other. For the lethal area com¬ 
putations, the total impact pattern of all classes of flechettes is divided 
into two parts: 
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• Areas in the ground plane which contain flechettes from only 
one class. 

• Areas in the ground plane which contain flechettes from two , 
or more classes. 

Lethal area computations are accomplished through a numerical 
integration of kill probabilities over all area. These kill probabilities 
can be easily calculated for areas in the grounci pattern which contain 
flechettes from only one class; however, those areas in which the pat¬ 
terns overlap present a more complex solution. The program initially 
assumes that no overlap occurs, and the area of each polygon (N) in 
all class patterns is individually calculated. The impact angle and 
terminal velocity of flechettes at each polygon vertex is calculated in 
the trajectory portion of the program. An average or mean value of 
velocity and impact angle is calculated for each polygon using these 
data. In addition, density of the flechettes is calculated by taking the 
total number of flechettes in the original polygon (determined'at release) 
and dividing by the polygon area at impact. Kill probability is calculated 
by: 

PK(N) = 1 - e PApPc (106) 


where P is flechette density in polygon N, Apis presented irea of the 
target, is conditional kill probability. A preliminary minition lethal 
area is calculated by summing the product of kill probability and area for 
all polygons: 


K 

AL = ^PpK(N) (A(N)) (107) 

1 

where K is the total number of polygons in all classes. 

At this point in the program a check is made to deter nine if the 
ground patterns from the various groups o '•■•lap. If only one class of 
flechettes is being run, or if no overlap oce-' . the munition lethal area 
is given by Equation (107). If two or more classes have flechettes falling 
in a common area, the overlap condition exists. Limits are then de¬ 
termined to localize the areas of overlap, and a rectangular grid of 
points is superimposed over this area. Each grid point is associated with 

68 







a prespecified area. The loss of lethal area (DALG) due to overlap is 
calculated for the incremental area associated with this grid point by: 



( 108 ) 


where ACRID is grid incremental area, PK(1,N) is kill probability in 
polygon (N) of class (I), and J is the total number of polygons which 
contain the grid point. 
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